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SUMMARY 

A unified model of voltage mode control (VMC) and current mode control (CMC) is proposed to 
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1 Introduction 



Occurrence of saddle-node bifurcation (SNB) [T] is generally unnoticed. A "different" instability and 
a possibility of two steady-state solutions in the buck converter under current mode control (CMC) 
in discontinuous conduction mode (DCM) have been reported [2l[3]. The "different" instability is 
actually a saddle-node bifurcation. In some specific instances when SNB is explicitly noticed [HE], 
no SNB boundary conditions have been derived. 

In this paper, a unified model of voltage mode control (VMC) and CMC is proposed to predict 
SNB. Based on the unified model, seven different approaches are applied to analyze SNB. Exact 
boundary conditions are derived based on these approaches. The boundary conditions define the 
critical values of converter parameters when SNB occurs. Applying different approaches actually 
leads to the same boundary conditions, further confirming the accuracy of the derived boundary 
conditions. The exact boundary conditions can be further simplified in various approximate closed 
forms for design purpose. 

The SNB may explain some sudden jumps or disappearances of steady-state solutions observed in 
DC-DC switching converters. Suppose SNB occurs at a critical parameter value. On one side of the 
critical parameter value (also the bifurcation point), there are multiple solutions (attractors) with 
different domains of attraction. When the system is perturbed, the state of the converter may jump 
from one attractor to another. On the other side of the critical parameter value, no solution exists. 
When the parameter crosses the critical value, the converters originally have multiple solutions, but 
now these solutions suddenly disappear. Those phenomena such as multiple solutions, sudden jumps 
and disappearances are undesirable in DC-DC switching converters. It is useful to understand the 
dynamics of SNB and its boundary condition to avoid its occurrence. 

The SNB actually exists in popular DC-DC converters, such as buck or boost converters, under 
various control schemes, such as VMC, CMC, or multi-loop state feedback control. The seven dif- 
ferent approaches are applied to accurately predict occurrence of SNB in buck and boost converters 
under these control schemes. Without loss generality, only continuous conduction mode (CCM) is 
considered. Analysis of SNB in DCM is reported separately [6]. 

The remainder of the paper is organized as follows. In Section [21 the operation of VMC and 
CMC is modeled in a single unified block diagram. In Section [3l steady-state analysis and small- 
signal analysis are presented. In Section UJ based on sampled-data slope-based analysis, general 
SNB boundary conditions for buck and boost converters are derived. In Sections 5-7, average, 
harmonic balance, and loop gain analyses are applied to predict SNB. In Sections 8 and 9, various 
approaches are applied to analyze SNB in buck and boost converters, respectively. Conclusions are 
collected in Section (THl 



2 Operation of Voltage/Current Mode Control 

The operation of a DC-DC switching converter under VMC or CMC can be described exactly by 
a unified block diagram model developed in (71 [8] shown in Fig. [H The control (reference) signal 
Vr controls the output voltage Vq in VMC. In CMC, Vr is generally denoted as ic to control the 
(peak) inductor current i/,. Denote the source voltage as Vs- In the model, Ai,A2 € R^^^, 
Bi,B2 £ R^^^, C,Ei,E2 G R^^^, and D G R'^^^ are constant matrices, where N is the system 
dimension. For example, N = 5 for a buck converter with a type III compensator. Within a clock 
period T, the dynamics is switched between two stages. Si and ^2. Switching occurs when the 
ramp signal h{t) intersects with the compensator output y := Cx -\- Du G R. Denote the ramp 
amplitude as V/i, and denote the switching frequency as fg ■= l/T and let ujs ■= 27r/s. 

Typical signal waveforms for VMC and CMC are shown in Fig. [2] and Fig. [3l respectively. In 
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Figure 1: Block diagram model for switching converter 
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Figure 2: Waveforms for voltage mode control 



Fig. [3l the ramp has positive slope, instead of negative slope as commonly seen in most literature, 
in order to be consistent with VMC. Other control schemes (average current mode control, for 
example) also fit the model of Fig. [TJ 



3 Steady-State and Small-Signal Analysis 
3.1 Periodic Orbits as Steady-State Solutions 

The periodic solution x^it) of the system in Fig. [1] corresponds to a fixed point x''(0) in the sampled- 
data dynamics. A typical periodic solution x^{t) is shown in Fig. HI where iP{d~) = Aix^ {d) + Biu 
and 3f{d^) = A2X^{d) + B2U denote the time derivative of x^{t) at t = d'^ and d"*", respectively. 
Let y^{t) = Cx^{t) + Du. In steady state, = CiP{t). Let the steady-state duty cycle be D 

and d := DT. Confusion of notations for capacitance C and duty cycle D with the matrices C and 
D can be avoided from the context. 



4 




I I I I Clock 

I Switch 

S-] S2 S-] S2 S-f S2 

Figure 3: Waveforms for current mode control 




Figure 4: A typical periodic solution x^{t) of a DC-DC converter in state space 



In steady state, 



Jo 



T-d 



From ([T]) and ([2]), one has 



(1) 
(2) 



x^{d) = (/-e^iV2(^-'^))-i(e^i'^ r \'^^"daB2U+ [^^''^daBiu) (3) 

Jo Jo 

Let Bi := [Bii,Bi2], B2 := [-B2i)-B22] to expand the matrices into two columns. The buck 
converter generally has Ai = A2 (invertible) , B21 = 0, and B12 = i?22- Then, from ([1]) and ([2]), one 
has 

x\d) = {I- e^i^)-!^^^^^"' - I)BiiVs - A^^Bi2Vr (4) 
The boost converter generally has Bi = B2, then 

f-T-d fd 

= (/-e^iV2(^-'^))-He^^W e^2"da+ / e^''' d(j)Biu := X{d)Biu (5) 

Jo Jo 
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3.2 Boundary Conditions Derived From Steady-State Analysis 

In steady state, 

Cx°((i) + Du- h{d) = (6) 

which is an equation in terms of d = DT. Generahy, if this equation has multiple solutions of d, 
SNB may occur if a converter parameter varies. When SNB occurs, only a single solution exists, 
which means that the curve Cx^{d) + Du — h{d) has a flat slope when SNB occurs. 

Take the buck converter for example, using ([4]), and differentiate ([6]) with respective to d, one 
has the SNB boundary condition, 

TC{I - e^'^r^e^^'^BuVs -Vh = (7) 
It will be shown later that applying different approaches also leads to the same condition. 

3.3 Small-Signal Analysis 

Using a hat " to denote small perturbations (e.g., Xn = Xn — x''(0)). From [71[8l[9], the linearized 
sampled-data dynamics is 

Xn+l = ^Xn (8) 

where the Jacobian matrix [H] 

yO{d-)-h{d) 

SNB occurs when one eigenvalue of $ is 1, and det[/ — $] =0. Since other bifurcations [10], such 
as pitchfork or transcritical bifurcations, also have one eigenvalue of $ at 1, the same analysis can 
be applied to analyze these bifurcations, omitted to save space. Although the general methodology 
of sampled-data analysis has been known in the last three decades, the closed form expression of 
([9]) was first derived and published in [3 [9], to the author's knowledge. All slope-based boundary 
conditions derived in this paper are based on the closed form expression of ([!]). 

4 Sampled-Data Slope-Based Analysis 
4.1 General Boundary Conditions 

Based on det[/ — <I>] = 0, the following theorem is obtained [3 p. 46]. 

Theorem 1 In a system as shown in Fig. 1, saddle-node bifurcation occurs when 

y\d-) - C{I - e-^2(^-'^)e-^i'^)-i(i;°(d-) - = h{d) (10) 

The proof is as follows. Suppose 1 is not an eigenvalue of QM{T-d)^A\d^ then 

det[/ - CD] = det[/ - e^2(T-d)gAid] ^^^^j ^ _ ^A2iT-d)^Mdyi^A2iT-d) i%d~) - ^°(^+) ^gAid] 

Cx^{d") — h{d) 

= det[I - e^2(T-d)g^i*ri ^ c^Md^j _ ^A2iT-d)^Mdyl^A2{T-d)A£)^Ad^^ 

yO{d-)-h{d) 

det [/ — $] = requires that the last term (inside of the second square brackets) of the last equation 
equals to zero, which leads to ([TO]) . 
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One can expand (jlOp in terms of x^{d), 

C{Aix\d) + Biu) - C{I - e-^2{T-d)g-Aid)-i(^^^ _ A2)x\d) + {Bi - B2)u) = h{d) (11) 

Note that the condition (|10p is valid for either VMC and CMC, and apphcable to general 
switching converters of any system dimension. Also note that in (llOp . the left side is in terms of 
the ripple slopes {y^{d~), x^{d~), and x^{d~^)), and the right side is the ramp slope h{d). As in the 
popular slope-based bomidary condition for subharmonic oscillation in CMC, the condition (jlOp is 
also slope-based. 

The condition (|10p can be proved to be equivalent to 

y\d+) - C(e^2(^-'^)e^i'^ - I)~\±\d") - ±\d+)) = h{d) (12) 

Using which one of (|10p or (jl2p depends on convenience. Since the proof for (jlOp is given, the proof 
for (fT2|) is omitted to save space. 



4.2 Buck Converter 

The buck converter generally has Ai = A2, B21 = 0, and B12 = -622- Using @, the boundary 
condition ([TT]) becomes 

C{I - e'^'^)-^e'^"^BuVs = h{d) (13) 
or in terms of -Ug, which shows the critical value of Vs when SNB occurs, 

*W (14) 



' C(/-e^i^)-ie^i'^Sii 

Both ([13]) and (fH|) are equivalent to d?]) which is based on the steady-state analysis. 

Based on the assumption that RC and y/LC are much larger than T, matrix approximations 
such as e^i^ I + AiT + AIt"^ /2 and {I + AiXy^ I-^iT can be applied. Then, the boundary 
condition (jl3p leads to 

1 Z 12 Vs 

or in terms of Vs, 

h{d) 



-^CA~'Bii + (i - D)CBii - (l^^q±Ml)CAiBnT 



(16) 



Remarks: 

(a) The left side of the boundary condition (|15p is a weighted combination of CA^^Bn/T, CBu 
and CAiBiiT. It will be shown that if the equivalent series resistance (ESR) Rc = 0, CAiBn 
dominates in VMC, while CBu dominates in CMC. For Rc > 0, either VMC or CMC has both 
the terms CBu and CAiBuT, indicating that the SNB conditions for VMC and CMC are closely 
related. 

(b) The boundary condition (jlSp seems to have a pattern. It is a hypothesis that the exact 
boundary condition has the form 

C( V 5n{D)A^T^)Bu = ^ (17) 
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where 5-i(-D) = -1, 5o{D) = (l-2D)/2, 5i{D) = {-1 + 6D - 6D'^)/12, etc., and d5n+i{D)/dD = 
6n{D). Further research on the series expression of (jl7p will be reported separately. However, since 
an exact condition as in (jl3p has been obtained, another exact condition in series expression may 
give additional insights but may be unnecessary. 

(c) The condition is exact, while (|15|) is approximate. Based on simulations, if the poles 
are smaller than w^/lO, the approximate condition (jlSp is close to the exact condition (jl3p . and it 
is generally adequate to predict SNB. For poles greater than Ws/10, using the exact condition (fT3]) 
gives more accurate results. 

(d) When the switching frequency is high (T is small), then the first term of psp dominates, 
and the boundary condition (jl5p becomes 

_ 1 CA-B„ » M (18) 

T Vs 

It will be shown later that this is the boundary condition derived from the average analysis. 

4.3 Boost Converter 

Let A((i) := I+{Ai-{I-e-^^^'^-'^'ie-'^''^)-^{Ai-A2))X{d) to simplify the equation. The analysis 
for the boost converter is similar to that for the buck converter. Using ([5]), the boundary condition 
(fTTj) becomes 

CA{d)Biu = h{d) (19) 

or in terms of Vg-, 

h{d) - CA{d)BuVr .... 

= — cmB^, — 

4.4 The "S plot": a Slope-Based Plot in the Real Domain 

Define an "S plot" as a function of D = d/T, 

S{D) := y\d~)-C{I -e-'^^^'^-'^h-'^"^Y^{±\d~)-±\d+)) (21) 
= C{I - e^i^)-^e^i'^e^i™t;s (for buck converter, from ^) (22) 
= CK{DT)Biu (for boost converter, from ([E])) (23) 

where, for large feedback gain, Vs ~ Vq/D for the buck converter and Vs ~ fo(l — D) for the boost 
converter. Note that (j2ip . (|22p . and (|23p are exact representations for the S plot, while (|15p is 
an approximate representation, which is generally accurate if the poles are smaller than w^/lO as 
discussed above. Then, from (jlOp . SNB occurs when 

S{D) = h{d) (24) 



5 State-Space Average Analysis 

Average analysis can be applied to analyze SNB when the switching condition is determined by the 
average signal values, such as in VMC, or in average CMC. If the switching condition is determined 
by the peak signal values, such as in peak CMC, the average analysis may be inaccurate. 

Let A = DAi + (1 — D)A2, and B = DBi + (1 — D)B2. In the state-space average model [3], 
the power stage dynamics is x = Ax + Bu, and the steady-state solution is 

- A-^'Bu := X (25) 
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Using a hat " to denote small perturbations, (e.g., x = x—X), the linearized open- loop dynamics 
[3] is 

£ = Ax + {{Ai - A2)X + {Bi - B2)u)D (26) 
In the closed loop, D = y/Vh = Cx/Vh- Then, the closed-loop dynamics is 

When SNB occurs, $a has an eigenvalue at 0, and det[<I>^] = 0. Suppose det[A] ^ then 

1 1 , r.n, r. A'^ ((Al - A2) X + (Bl - B2)u)C , , , 

det[#A] = det[yl] det[/ + ^' ^ i±A_-\ = q (28) 

leading to the following boundary condition 

Vh + CA~^{{Ai - A2)X + {Bi - B2)u) = (29) 

The buck converter generally has Ai = A2 (invertible) , B21 == Oatxi, and B12 = B22, and 
the boundary condition ([29|) becomes (equivalent to (flSj) based on the sampled-data slope-based 
analysis) 

Vh + CA-^BuVs = (30) 
The boost converter generally has Bi = B2, and the boundary condition ()29p becomes 

Vh + CA~\Ai-A2)X = (31) 

Note that the condition (I30p is independent of D, and is only an indication of the closeness to SNB, 
while the condition (I3ip is a function of D through A = DAi + (1 — D)A2- 

6 Harmonic Balance Analysis 

Applying harmonic balance (HE) analysis to predict subharmonic oscillation has been reported in 
[71 [11]. Here, HE analysis is applied to predict SNB. Consider a buck converter power stage with a 
control-to-output (D-to-Vo) transfer function G^dis)- In the converter, there is an ON switch and 
an OFF switch (sometimes substituted by a diode). Let the voltage across the OFF switch (or 
the diode) be Vd- The waveform of Vdit) is a square wave with the high voltage at Vg and the low 
voltage at 0, which can be represented by Fourier series (harmonics) 

00 

Vd{t)= V c„e^""^* where c„ = ^(1 - e"^"'^^'^) (32) 



n=—oo 



The HE analysis can be applied to VMC, or CMC with voltage loop open/closed. 
6.1 VMC 

Let the Vd-to-Vo transfer function be Gv{s). One has [31 p. 470] 

G.is) = ^ = ^RcC + l 
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Let the compensator transfer function (from Vq to —y (negative sign due to negative feedback)) 
be Gc{s). Let the transfer function from to —y be G{s) = Gcis)Gy{s) = Gc{s)Gi,d{s)/vs- For 
designation purpose, G{s) is called a Harmonic Balance (HB) gain, which is proportional to the 
loop gain T{s) = Gc{s)Gyd{s)/Vh by 

G{s) = y^T{s) (34) 

Vs 

Then, the signal y{t) at the output of the compensator is 

y{t) = Vr+Gc* {vr - Vo) (35) 

= {l + Gcmvr-Gc*Vo (36) 

= {l + G,{0))vr-GcG,*Vd{t) (37) 

= {1 + Gcmvr-G^vd{t) (38) 

oo 

= {1 + G,{0))vr- Cne^'^'^^'GijnoOs) (39) 

n=— oo 

oo 

= (1 + Gc{0))vr - VsDG{0) - 2Re ^ c„e^"'^»*G(jna;,) (40) 

n=l 

where Re denotes taking the real part of a complex number, and -k denotes convolution. 

The intersection of y{t) and h{t) determines the duty cycle and hence the waveform of Vci{t). 
By "balancing" the equation y[t) = h{t) (written in Fourier series form) at the switching instants, 
conditions for existence of periodic solutions and SNB can be derived. 

In steady state, 

oo 

y{d)-h{d) = {l + Gcmyr- Cne^'"'"'^G{jnUs)-h{d) = Q (41) 

n=— oo 

which is an equation in terms of d = DT. When SNB occurs, only a single solution exists, which 
means that the curve y{d) — h{d) has a flat slope at the critical value of d when SNB occurs. 
Differentiate (j4ip with respective to d, one has the following theorem. 

Theorem 2 Given a closed-loop buck converter with a control-to- output transfer function Gyd{s) 
and a compensator transfer function Gc{s), let the HB gain G{s) = Gc{s)Gyd{s) /vg- SNB occurs 
when 



oo 
n=—oo 



Y e^^^^''GijnoJs) + Vh = (42) 



The condition (j42|) can be proved to be equivalent to d?]) based on the steady-state analysis or 
(jlSP based on the sampled-data slope-based analysis, showing that applying different approaches 
leads to the same condition. The boundary condition (j42p can be expressed in various forms for 
design purpose. For example, (|42p leads to any of the following conditions, 

oo 

VsG{0) + 2vsReY(^^^''"^G{jnUs) + Vh = (43) 

n=l 

Re[f;e^-2-^G(jn^,)] = -YH^^^ (44) 



2vs 

n=l 
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or in terms of Vs, 

= G(0) + 2ReE~ie^'2"-^G(inu;,)] ^^^^ 

(47) 

Generally, Gy{s), Gc{s) and thus G{s) are low-pass filters. The summation in (j44p can be 
approximated by the term that involves G{s) with the smallest argument. Therefore, (j44p becomes 

Re[e^--G(,a..)] = -^^^^ (48) 
The "H plot": a Nyquist-like plot in the complex plane. 

Note that the left side of (I44p is a function of D, oog, and the HB gain G{s), where G{s) is further 
a function of the power stage and compensator parameters. Let 

oo 

if(I)):=5^e^2"-^G(ina;,) (49) 



n=l 



Then, SNB occurs when 



Re[i^(«)l = .}^±^ ,50) 



For designation purpose, H{D) is called an H plot because it is similar to the Nyquist plot in the 
complex plane for design purpose. Given a desired range of D, one can plot H{D) according to 
(I49p to determine whether SNB occurs in this range of D. 



6.2 CMC with Voltage Loop Open 

In CMC with the voltage loop open, the controlled output is the inductor current ii. Let the peak 
inductor current control signal be v (also denoted as Vr in Fig. [T]) which controls the peak inductor 
current. One has y{t) = ic{t) — iL(t). Similar to (j33|) . the Vd-to-ii transfer function is [Sj p. 470] 

G (s) ■= + + ^ (51) 

LC(1 + f )s2 + (^ + ii,C)s + 1 ^ ' 

Since no extra compensator (except the compensating ramp h{t)) is added in the current loop, 
Gc{s) = 1, and one has the HB gain G{s) = Gc{s)Gi{s) = Gi{s). Since the expression of G{s) is 
derived, the rest of analysis is similar as in VMC. 

6.3 CMC with Voltage Loop Closed 

Here, there are two feedback loops. Let Gc{s) be the voltage loop compensator. One has y{t) = 
ic{t) — iiit) = Gc* {vr — Vo{t)) — iL{t) = Gc{Q)vr — {Gc * G^, + Gi) -k Vd{t). Then, the HB gain is 

G{s) = G,{s)G,{s) + G,{s) (52) 

Since the expression of G{s) is derived, the rest of analysis is similar as in VMC. 
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7 Prediction of SNB Based on Loop Gain 

Since the HB gain G{s) is proportional to the loop gain T{s) = G{s)vs/Vh, (02]) directly leads to 
the following theorem 

Theorem 3 Given a closed-loop buck converter with a loop gain T{s), SNB occurs when 

J2 e^'2"-^r(jn6.,) = -1 (53) 

Tl= — oo 

The condition ()53p can be expressed in various forms, 

oo 

T(0) + 2Re[^e^'2"^^r(jna;,)] = -1 (54) 



n=l 

oo 



Re[^e^-2-^r(jn...)] = -li^ (55) 

n=l 

The summation in (j55p can be approximated by the term that involves T(s) with the smallest 
argument. Therefore, (|44p becomes 

Re[e^'-^T{jUs)] - -I^^ (56) 
It should be noted that (|56p is only an approximate condition, and the exact condition is ()53p . 
The "LI plot" in the real domain. 

Note that the boundary condition (I55p is a function of D, ojs, and the loop gain T(s), where T{s) 
is further a function of Vg, Vh, tlie power stage and compensator parameters. Define an LI plot, 
which is a real function, as 

oo 

L,{D) := e^^'^^^'TUkoOs) (57) 

fc=— 00 

Then, SNB occurs when 

Li{D) = -1 (58) 

The "L2 plot" in the real domain. 

Since in some situations, = 0, such as in CMC with no ramp compensation, T(s) = G{s)vs/Vh 
becomes infinite. In such a situation, define an L2 plot, which is also a real function, as 

00 

L2{D):= e^'"-''G(jn^,) (59) 

n=— 00 

Then, from (j42p . SNB occurs when 

L2{D) = (60) 

Vs 

Since (j58p . (|60p and (jl3p are exact conditions for occurrence of SNB, they are equivalent, and the 
LI and L2 plots have matrix forms 

L,{D) = = _rC(/ - e^^'^r^e^'^Bn (61) 
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8 Buck Converter 



This section shows occurrence or non-occurrence of SNB in the buck converter under various control 
schemes. Boundary conditions for specific schemes are derived if SNB occurs. Similar analysis is 
applied to the boost converter in the next section. 

8.1 CMC with Voltage Loop Open 

Let the state x = {iLj^c)', where is the inductor current, and vc is the capacitor voltage. Then, 



Bi 
El 



A. 



? t 

C RC 



[-611,-812] = 

E2 = [Q 1 



I 




Bo 



[B21 , B' 



22 








(62) 



In CMC with the voltage loop open, a switching occurs when — iiit) = h{t). One has 
y = ic — and 

C=[-1,0] I? =[0,1] 

CA~^Bii = ^, CBu = ^ CAiBu = 

By simple algebra, the boundary condition ([15]) leads to 



j,_l + K ^ Lh{d) 



(63) 



where K := 21/ RT as defined in [3]. Since D <l, then from ([63]), SNB in the CMC buck converter 
occurs only when K <\. Also from ()63p . SNB occurs when the converter is unstable {D > 1/2). 

The boundary condition (j63l) can be also derived by the harmonic balance analysis. Based on 
the facts that for < D < 1 



^sin( 



k=l 



2TTkD^ 
k ' 



and using (|5T|) . the boundary condition (p3|) becomes 

00 



n=l 



R 



vJT 
L 



fc=i 



,2-KkD, , , 

sin(— -) + Vh (64) 



which leads to 
condition. 



(65) 

It is interesting to note that completely different approaches lead to the same 



Example 1. Consider a CMC buck converter with the following parameters: fs = 5 V, /s = 200 
kHz, i? = 5 ri, L = 5 /iH, and C = 40 ^F. Here, K = 0.4 and the converter is at light loading and 
has both CCM and DCM solutions. 

Let ic be the bifurcation parameter. The bifurcation diagram is shown in Fig. [5j SNB occurs 
when ic = 1.225, D = 0.7, and Vq = 3.5. Note that the focus of this paper is on the T-periodic 
solution. Those non- T-periodic attractors are not shown in the bifurcation diagram to prevent 
detraction of the focus on SNB. 
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Figure 5: Bifurcation diagram showing SNB at ic = 1.225 and CCM-DCM transition at ic = 1-2, 
solid line for stable DCM solution, dotted line for unstable CCM solution 
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Figure 6: Coexistence of two unstable T-periodic CCM solutions (dashed line) and a stable 2T- 
periodic solution (solid line), ic = 1-21 



For ic > 1.225, there is no solution. For v G (1-2, 1.225), there coexist at least three solutions: 
two unstable (with D > 1/2) CCM solutions and a stable attractor. Take v = 1-21, for example. 
The two unstable CCM solutions, with duty cycles 0.62 and 0.78, and the stable attractor are 
shown in Fig. [H The stable attractor is 2T-periodic, in CCM for the first period and in DCM 
(with a tiny duration when the inductor current is zero) for the second period. Note that the three 
solutions have the same peak inductor current. 

As ic is increased to 1.223, the two unstable CCM solutions move closer, with duty cycles 0.67 
and 0.73, as shown in Fig. [71 At v = 1.225, SNB occurs when the two unstable CCM solutions 
merge together. 

At ic = 1.2, the converter transitions from unstable CCM to stable DCM through a CCM-DCM 
border-collision bifurcation [12]. For ic G (1, 1-2), there exist one stable (with conversion ratio [3] 
M < 2/3) DCM and one unstable CCM (with D > 1/2) solutions. For ic < 1, there exists only 
one stable DCM solution. 

The SNB boundary conditions will be derived using various approaches. Here, K = 0.4 and 



14 



2 3 
Time (Second) 



4 5 
X 10"° 



3.3 3.35 3.4 3.45 3.5 3.55 3.6 3.65 3.7 
Output voltage 



(a) (b) 
Figure 7: Coexistence of two unstable CCM solutions, ic = 1.223, (a) time domain, (b) state space 
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Figure 8: The S plot, S{D) 



h{d) = 0, then from ([63]) . SNB occurs at D = 0.7. The S plot (based on the slope-based analysis), 
H plot (based on the harmonic balance analysis), L2 (based on the loop gain analysis) are shown 
in Fig. [8l Fig. [9l and Fig. [TOl respectively. They all show that SNB occurs exactly at D = 0.7. □ 



8.2 VMC: No SNB 

Without loss of generality, let the voltage loop has a proportional feedback gain kp. One has 
y = kp{vr — Vo) and 

C=[{),-kp\ D=[{),kp\ 
CA-^Bn = kp, CBu = CAiBu = ^ 

The boundary condition p3|) becomes 

T\1-6D + 6D^) Vn 
+ V2LC =k^s ^^^^ 
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Figure 9: The H plot, H{D) 
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Figure 10: The L2 plot, L2{D) 



Generally, < 12LC, and the left side of is negative while the right side is positive. The 
boundary condition is not met, and SNB does not occur in this situation. 

8.3 Multi-Loop State Feedback 

Consider a buck converter under multi-loop state feedback, y = Vr — kii^ — k^vc- One has 

C = -[ki,k^] -D=[0,1] 
CA-^Bn = I + K„ CBn = ^ C^i^ii = ^ 

The boundary condition (jlSp leads to 



^^(z,-iI±i) + ^^(-i + ^!(^5|±^,.AM ,67) 

where the left side of (j67p is also an approximate form for the S plot. For <^ 12LC, solving (|67p 
for D, SNB occurs when 

K + 1 Lh(d) Lky , , 
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Figure 11: Bifurcation diagram showing stable (solid) and unstable (dashed) solutions. SNB occurs 
at Vs = 20 and D = 0.7. 



Example 2. Consider a buck converter under multi-loop state feedback [H p. 232]. In [T], a digital 
control is used, but can be approximated as a multi-loop analog control, with y{t) = Vr — kiii — kyVc 
intersecting with h{t) to determine the duty cycle. The converter parameters are T = 400 /is, L = 20 
mH, C = 47 fiF, R = 22 Q, ki = 2.1435, k^ = -0.1383, Vr = 0.2152 V, and Vh = 1. 

Using Vs as the bifurcation parameter, the bifurcation diagram is shown in Fig. [TT\ showing 
that SNB occurs at Vs = 20 and D = 0.7. Generally in SNB, there is a hysteretic loop as shown 
in Fig. [TT] In the figure, the upper solid line is for the operation when the switch is always on 
(hence D = 1), and the dashed line and the lower solid line are for unstable and stable T-periodic 
solutions respectively with duty cycle less than 1. For Vs between 19.25 V and 20 V, there are three 
solutions: one stable T-periodic solution, one unstable T-periodic solution, and the third (stable) 
solution being that the switch is always on. When the converter operates with a periodic solution 
and Vs is increased a little from 20 V, the output voltage will jump up from 14 V to 20 V. Similarly, 
when the converter operates with D = 1 and Vs is decreased a little from 19.25 V, the output 
voltage will jump down from 19.25 V to 10 V. The jumping up and down forms a hysteretic loop. 
The S plot and its approximate closed form (I67p are shown in Fig. \T2\ both showing that SNB 
occurs exactly at T* = 0.7. Also, from ([681) . using the fact that SNB occurs at Vs = 20, then ([68]) 
gives D = 0.71, which is close to the exact critical value. □ 



9 Boost Converter 
9.1 VMC 

Consider a boost converter under VMC with proportional feedback gain kp. One has y = kp[vj.—vc)- 
Let the parasitic resistance associated the inductor be r. Then, 
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Figure 13: The curve of vc{D)/vs, for p = (dashed hne) and p = 0.05 (solid line). 



Since the average of y is used in VMC to determine the duty cycle, the average analysis can be 
applied. Let p = r/R and k = kp/Vh- From (p5]) . 



X = -A-^Bu 



2_ 

R 

1-D 



vc{D) 



(70) 



The curve of vc{D)/vs is shown in Fig. 1131 for p = Q and p = 0.05. For p = 0, the curve of vc{D)/vs 
increases monotonously as D increases, and one output voltage has only one corresponding duty 
cycle. For p = 0.05, the curve of vc{D)/vs is A-shaped [3l p. 45]. One output voltage has two 
corresponding duty cycles. The SNB is generally related to coexistence of multiple solutions [T]. It 
will be shown that SNB occurs only when p > 0. 

Using ()69p . the boundary condition ()3ip based on the average analysis leads to 



(p + (1 - Dyy + KVsiii - DY -p)) = Q 



(71) 
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Figure 14: Bifurcation diagram showing stable (solid) and unstable (dotted) solutions. Hopf bifur- 
cation occurs at Vr = 4.92, and SNB occurs at Vr = 7.1. 



Solving ()7ip for D, then SNB occurs when 



If the parasitic inductor resistance r = (hence p = 0), then SNB does not occur for D <1. This 
shows a dramatic difference if the parasitic inductor resistance r is not considered in the model. 
If the parasitic resistance is not modeled, one may be misled by the wrong dynamics (no SNB) 
and steady-state solutions (only one solution), while the actual circuit (with r > 0) has SNB and 
multiple solutions. 

Example 3. Consider a boost converter under VMC with kp = 2. The converter parameters 
are = 3 V, = 1, fs = 600 kHz, L = 1 /iH, C = 100 /iF, R = 2 Vt, and parasitic inductor 
resistance r = 0.1 0. One has p = r / R = 0.05. The bifurcation diagram is shown in Fig. [TH The 
Hopf bifurcation [13] occurs at Vr = 4.92 and SNB occurs at Vr = 7.1. 

Note that the focus of this paper is on the T-periodic solutions. Let Vr = 7, for example. The 
two unstable T-periodic solutions are shown in Fig. [T5l One solution has D = 0.74, and the other 
has D = 0.81, agreed with Fig. [Ml Other stable attractors may coexist with the two unstable 
solutions. In this particular example, the converter is actually bistable for 4 < t;^ < 4.92 and 
monostable for Vj. > 4.92. The other stable attractor is {iL,vc) = (vs/rjO) = (30,0), and D = I 
(one switch is always on) because y{t) = kp{vr — vc) = kpVr > 1 > h{t). Such an attractor does 
not exist if r = because it requires i_L = Vs/r = oo. With different converter parameters, the 
attractor may be quasi-periodic [T], associated with the Hopf bifurcation [13]. Other attractors 
associated with the border-collision bifurcation [H] may also exist when y(t) of the attractor is out 
of the bounds of h{t). These non- T-periodic attractors are not shown in the bifurcation diagram 
to prevent detraction of the focus on SNB. 

From ()72p . SNB occurs when D = 0.78, the prediction of the SNB agrees with the bifurcation 
diagram. The S plot (Fig. I16|) shows that SNB occurs exactly at D = 0.78. □ 



D 




(72) 
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-y/p (for large k) 



(73) 
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Figure 15: Ramp h{t) (solid line) and two unstable T-periodic solutions y^(t) (dashed line), coex- 
isting with stable attractor {iL,vc) = (30,0) 
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Figure 16: The S plot, S{D) 



9.2 CMC with Voltage Loop Closed 

In CMC with the voltage loop closed, the voltage loop output controls the peak inductor current. 
Without loss of generality, let the voltage loop has a proportional feedback gain kp. In steady state, 

ic-iL = kp{vr - Vo) -iL = h{t) (74) 

The model is the same as in ([69]) . except that here C = [—1,—kp]. 

As in VMC, one can prove that SNB occurs when p > 0, and SNB does not occur when p = 0. 
In steady state, (fTll) is rearranged as 

Vr=Vo^ T (75) 

For large /cp, (i_L + h[d))/kp can be ignored. Since the curve Vq as a function of D is A-shaped as 
shown in Fig. [13] if p > 0, one value of Vr corresponds to two duty cycles and SNB occurs. 

In contrast, if p = 0, fo increases monotonously as D increases as shown in Fig. [131 one value 
of Vr corresponds to only one duty cycle and SNB does not occur. 
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Figure 17: Bifurcation diagram showing stable (solid) and unstable (dotted) solutions. SNB occurs 
at Vr = 17.71 and D = 0.91, and period-doubling bifurcation occurs at Vr = 8.2. 




Figure 18: The S plot, S{D) 



Example 4. Consider a boost converter under CMC with no ramp compensation (Vh = 0). The 
voltage loop has a feedback gain kp = 2. The converter parameters are the same as in Example 3. 

The bifurcation diagram is shown in Fig. \T7\ showing that SNB occurs at Vr = 17.71 and 
D = 0.91, and the period-doubling bifurcation [13] occurs around D = 0.5. The S plot (Fig. [T8]) 
shows that SNB occurs exactly at D = 0.91. □ 



9.3 CMC with Voltage Loop Open: No SNB 

In CMC with the voltage loop open, the current control signal ic controls the peak inductor current. 
From (|69p . the inductor current ripple AIl = (vs — rlL)DT / L. The peak inductor current Ipeak = 
II + increases monotonously as D increases. Therefore, given a value of current control 

signal ic, only one solution exists and SNB does not occur. 
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Figure 19: Bifurcation diagram showing stable (solid) and unstable (dashed) solutions. Saddle-node 
bifurcation occurs at Vr = 0.496 and D = 0.65. 



9.4 Multi-Loop State Feedback 

Consider a boost converter under multi-loop state feedback, y = Vr — kflL — kyVc- The model is 
the same as in (j69p . except that here C = —[ki,ky] and D = [0,1]. The sampled-data analysis is 
similar to the case for the boost converter under VMC, omitted here to save space. The average 
analysis can be also applied to obtain an approximate boundary condition. From (I3ip . the boundary 
condition leads to 

Since D < 1, SNB does not occur if both ki and k^ are positive because the condition (j76|) is not 
met. SNB occurs only when at least one of k^ and k^ is negative. 

Example 5. Consider a boost converter under multi-loop state feedback [15\ p. 90]. The 
converter parameters are T = 2 /xs, = 4 V, L = 5.24 /iH, C = 0.2 fiF, R = 16 i}, ki = —0.1, 
ky = 0.01, Vr = 0.48, and Vh = 1. This converter has been analyzed with different views. In |15j . 
this converter is used to illustrate the inaccuracy of the average model about the converter stability. 
It is later proved in [71 p. 75] that the average model is actually accurate about local stability, but 
the converter is not globally stable. It is then shown in [16] that two solutions (one stable and the 
other unstable) coexist, but without linking multiple solutions with SNB. 

The bifurcation diagram is shown in Fig. [T9l showing that SNB occurs at Vr = 0.496 and 
D = 0.65. For Vr = 0.48, there are actually three solutions: one stable T-periodic solution, one 
unstable T-periodic solution, and the third (stable) solution being that one switch is always on 
(which does not exist in the real circuit because infinite inductor current is required). The S plot 
(Fig. [20]) shows that SNB occurs exactly at D = 0.65. Applying the average analysis can also give 
approximate results. Solving (j76]) gives D = 0.685, which is close to the exact critical value. □ 

10 Conclusion 

A unified model of VMC/CMC is proposed to predict SNB. In the unified model, confirmed with 
simulations and the exact sampled-data analysis [TJI^, SNB boundary conditions for VMC/CMC 
are derived and they have similar forms for VMC/CMC. The boundary conditions are expressed 
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Figure 20: The S plot, S{D) 



in terms of signal slopes, or other converter parameters. The sloped-based boundary conditions 
derived are exact, and can be further simplified in various approximate closed forms for design 
purpose. The sloped-based boundary conditions are applied to analyze buck and boost converters 
under various VMC/CMC control schemes. Since all the boundary conditions are derived in closed 
forms, one can see the effects of various converter parameters on the occurrence of SNB. 
Many approaches can be applied to analyze SNB, including 

1. simulation to plot the bifurcation diagram 

2. average analysis 

3. sampled-data small-signal analysis using the Jacobian matrix to determine the pole location 

4. steady-state analysis to determine whether multiple solutions exist 

5. sampled-data slope-based analysis and using the S plot 

6. harmonic balance analysis and using the H plot 

7. loop gain analysis and using the L1/L2 plot 

The first three approaches are known, but no closed form boundary conditions have been derived in 
the past research. The last four approaches are first proposed in this paper to derive the boundary 
conditions to the author's knowledge. Applying different approaches lead to the same boundary 
conditions, providing convincing evidences about the accuracy of the derived conditions. Each 
approach has its own merits, and complement the other approaches. For example, the slope- 
based analysis analyzes the converter in the time domain and expresses the boundary conditions in 
matrix forms. The harmonic balance analysis analyzes the converter in the frequency domain and 
expresses the boundary conditions in terms of harmonics of the switching frequency. Since most 
power stage and compensator transfer functions are low-pass filters, the harmonic balance analysis 
is particularly useful to analyze the converter with poles close to the switching frequency. 

Similar to the Nyquist plot which predicts stability for the continuous-time system, many new 
plots are proposed to predict SNB. The S plot is in matrix form, while the other plots are in terms 
of harmonics. The H plot is in the complex plane, while the other plots are in the real domain. All 
of these plots are exact and can accurately predict SNB. These plots also have simplified forms for 
design purpose. 
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Coexistence of multiple steady-state solutions is generally related to SNB, and steady-state anal- 
ysis can predict SNB. Generally, SNB occurs when two solutions (with different duty cycles) coexist 
corresponding to the same value of a bifurcation parameter. Therefore, the shape (monotonously 
increasing or A-shaped) of a controlled signal in steady state as a function of duty cycle deter- 
mines whether SNB occurs. For example, in the CMC buck converter, the peak inductor current 
is controlled and its steady-state curve as a function of duty cycle is A-shaped. Two steady-state 
solutions of different duty cycles have the same value of peak inductor current, and SNB occurs in 
the CMC buck converter. For the boost converter, the steady-state output voltage as a function 
of duty cycle is A-shaped (see Fig. [T1I|) if the parasitic inductor resistance is modeled. Therefore, 
SNB occurs in the boost converter if the output voltage is controlled, such as in VMC, or in CMC 
with the voltage loop closed. However, if the parasitic inductor resistance is zero, the steady-state 
output voltage monotonously increases as the duty cycle increases (see Fig. [T3|) . SNB does not 
occur in VMC, or in CMC with the voltage loop closed if the parasitic inductor resistance is zero. 
Similarly, in the boost converter, the steady-state inductor current monotonously increases as the 
duty cycle increases, and SNB does not occur if the inductor current is controlled, such as in CMC 
with the voltage loop open. 
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